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On  Formulating  a  Simple  D-Region  Model 


1.  INTRODUCTION 

The  D-region  is  the  smallest  portion  of  the  ionosphere  by  far.  The  vertical 
electron  content  of  this  layer  is  less  than  0.  01  percent  that  of  the  entire  ionosphere, 
which  is  nominally  6X  lO^^  cm'^  at  mid-latitudes  during  daytime.  Even  for  very 
disturbed  conditions,  the  fractional  content  of  the  D-region  does  not  exceed 
1.  0  percent.  Nevertheless,  the  D-region  is  responsible  for  LF  and  VLF  propaga¬ 
tion,  "survivable"  communication  modes.  Furthermore,  radio  wave  absorption 
can  be  appreciable  in  the  disturbed  D-region  at  higher  frequencies  like  HF  and  VHF. 
It  is  well  known  for  example,  that  HF  radio  waves  can  be  totally  absorbed  during  a 
Solar  Proton  Event  (SPE),  hence,  the  earlier  common  name.  Polar  Cap  Absorption 
(PCA)  event. 

The  large  neutral  but  low  charged  particle  concentrations  of  the  D-region  have 
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made  observations  of  this  region  difficult.  Thrane  and  Sechrist  have  discussed 
these  difficulties,  plus  the  uncertainties  associated  with  the  various  types  of 
electron  concentration  measurements. 


(Received  for  publication  20  October  1987) 

1.  Thrane,  E.  V.  (1974)  Methods  of  measurements  and  results  of  lower  ionospheric 

structure,  Akademie-Verlag,  Berlin,  pp.  3-21. 

2.  Sechrist,  C.F.,  Jr,  (1974)  Comparison  of  techniques  for  measurements  of 

D-region  electron  densities.  Radio  Sci.  137- 149. 
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The  D-region,  with  its  attendant  large  neutral  gas  concentrations  is  much 
more  difficult  to  model  than  the  other  (higher)  layers  of  the  ionosphere.  Three- 
body  processes  are  important.  The  composition  of  negative  and  positive  ions  is 
not  simply  related  to  the  major  gases,  because  minor  species  can  play  a  role  in 
three-body  ion  chemistry.  (The  chemistry  of  the  higher  ionospheric  layers  is 
limited  to  binary  processes. ) 

Models  of  the  D-region  prior  to  about  1965  dealt  with  a  very  limited  number 
of  ionic  species,  like  O^,  NO^  and  ©2.  The  pioneering  positive  ion  measurements 
of  the  D-region  [Narcisi  and  Bailey^]  showed  that  the  dominant  species  were  proton 
hydrate  (PH)  ions,  H„0"'’  *  (H„0)  .  Major  negative  ions  include  CO~,  HCOl  and 

NO,  [Arnold  et  al;  Narcisi  et  al^].  The  degree  of  hydration  of  these  ions  in  the 

^  6 
D-region  is  not  great,  especially  for  the  first  two  of  these  ions  [Keesee  et  al  ]. 

Temperature  and  water  vapor  concentrations  are  key  factors  in  the  hydration  level 

of  these  ions. 

Models  of  the  D-region  are  very  complex  since  the  distribution  of  minor 

species  like  NO  and  H2O  must  be  known.  Temperatures  are  also  required.  Many 

of  the  ionic  reactions  depend  strongly  upon  the  temperature.  Other  minor  species 

that  participate  in  the  chemistry  of  the  D-region  include  O,  O^*  NO2,  and  CO2. 

Global  models  of  these  species  and  the  temperature  field  are  extremely  difficult 
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to  formulate,  although  Brasseur  and  DeBaets  have  constructed  a  2-D  model  of 
the  D-region.  However,  an  accurate,  predictive  model  of  the  D-region  is  difficult 
to  formulate  and  is  cumbersome  since  so  many  parameters  are  required.  Constant 
monitoring  of  the  various  minor  species  is  impractical,  although  real  time  monitor¬ 
ing  of  energetic  precipitating  particles  and/or  x-rays  must  be  undertaken  for  the 
modelling  of  disturbed  events. 

3.  Narcisi,  R.  S. ,  and  Bailey,  A.D.  (1965)  Mass  spectrometric  measurements  of 

positive  ions  at  altitude  from  64  to  112  kilometers,  J.  Geophys.  Res. 
20:3687-3700.  - 

4.  Arnold,  F, ,  Kissel,  J.  ,  Krankowsky,  D. ,  Wieder,  H. ,  and  Zahringer,  J. 

(1971)  Negative  ions  in  the  lower  ionosphere;  A  mass  spectrometric 
measurement,  J.  Atmos.  Terr.  Phys.  33:1169-1175. 

5.  Narcisi,  R.  S. ,  Bailey,  A.  D. ,  Wlodyka,  L.  E.  ,  and  Philbrick,  C.R.  (1972) 

Ion  composition  measurements  in  the  lower  ionosphere  during  the  November 
1966  and  March  1970  solar  eclipses,  J.  Atmos.  Terr.  Phys.  34:647 -668. 

6.  Keesee,  R.  G. ,  Lee,  N. ,  and  Castleman,  A.W.,  Jr.  (1979)  Atmospheric 

negative  ion  hydration  derived  from  laboratory  results  and  comparison  to 
rocket-bo  me  measurements  in  the  lower  ionosphere,  J.  Geophys.  Res. 
24:3719-3722.  - - 

7.  Brasseur,  G.  ,  and  DeBaets,  P.  (1986)  Ions  in  the  mesosphere  and  lower 

thermosphere:  A  two-dimensional  model,  J.  Geophys.  Res. 

91:4025-4046. 


Based  upon  all  the  knowledge  gained  over  the  last  two  decades  particularly, 
and  in  view  of  the  complexities  for  physically  modelling  the  D-region,  it  will  be 
argued  here  that  this  region  may  be  equally  well  and  much  more  simply  modelled 
by  the  use  of  a  bulk  electron  loss  coefficient  for  the  daytime  and  an  analytical 
expression  for  nighttime.  This  approach  is  probably  best  for  predictive  purposes. 
This  argument  derives  from  a  considerably  rapid  ion  chemistry  and  a  low  level  of 
ionization,  that  is,  a  fairly  "frozen"  chemistry  [Swider  ].  Charged  particle  con- 
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centrations  never  approach  more  than  about  10  the  total  neutral  particle 
concentrations . 


2.  EFFECTIVE  RECOMBINATION  COEFnCIENT  OF 
THE  DAYTIME  D-REGION 

The  effective  recombination  coefficient  of  the  D-region  as  derived  from  the 
continuity  equation  for  electrons  under  steady  state  conditions  and  in  the  absence 
of  transport  processes  is; 


4*  =  (Qq  +  Ao.)(l+X)  =  q'[e]‘ 


(1) 


where  ajj  is  the  mean  dissociative  recombination  coefficient,  o.  is  the  mean 
ion-ion  recombination  coefficient,  X  is  the  ratio  of  negative  ions  to  electrons,  q  is 
the  total  ion-electron  pair  production  rate,  and  [e]  is  the  electron  concentration. 
The  notation  of  Adams  and  Masley,  SI?  for  has  been  adopted.  They  suggested 

that  ^  varies  only  with  altitude.  It  appears  that  this  suggestion  is  most  likely  to 
be  true  if  minor  neutral  species  are  not  involved  with  the  ion  chemistry,  that  is, 
if  all  chemistry  is  frozen  (fast).  However,  does  depend  somewhat  on  q,  [H  2O], 
and  temperature.  (Some  changes  in  certain  significant  minor  neutral  species 
concentrations  are  believed  to  occur  for  large  ionization  production  rates  [Swider 
and  Keneshea,  Solomon  et  al^  ^]. ) 


8.  Swider,  W.  (1988)  Electron  loss  and  the  determination  of  electron  concentra¬ 

tions  in  the  D-region,  PAGEOPH,  in  press. 

9.  Adams,  G.  W.  ,  and  Masley,  A.J.  (1965)  Production  rates  and  electron  densities 

in  the  lower  ionosphere  due  to  solar  cosmic  rays,  J.  Atmos.  Terr.  Phys. 
^7:289-298.  - 

10.  Swider,  W,  ,  and  Keneshea,  T.J.  (1973)  Decrease  of  ozone  and  atomic  oxygen 

in  the  lower  mesosphere  during  a  PCA  event.  Planet.  Space  Sci. 
^1:1969-1973.  - 

11.  Solomon,  S.  ,  Reid,  G.C.,  Rusch,  D.  W. ,  and  Thomas,  R.  J.  (1983) 

Mesospheric  ozone  depletion  during  the  solar  proton  event  of  July  13,  1982, 
Part  II.  comparison  between  theory  and  experiment,  Geophys.  Res.  Lett. 
10:257-260.  - 


For  the  quiet  D-region,  the  photo-ionization  of  NO  by  HLya  is  the  main  source 
of  electrons  in  the  daytime.  A  rapid  sequence  of  temperature-dependent  three- 
body  and  "switching"  reactions  leads  from  NO^  ions  to  NO"*^*  (H20)2  ions  and  thence 
to  proton -hydrate  (PH)  ions,  (H20)^.  with  the  n  =  2  ion  formed  first.  Thus, 

water  concentrations  must  be  known  as  well  as  temperatures. 

For  disturbed  conditions,  ot  ions  are  the  principal  precursor  ions  for  PH 
ions.  (The  more  abundantly  formed  N2  ions  are  unobservable  since  charge 
transfer  with  0„  is  virtually  instantaneous  at  D-region  altitudes.  )  The  sequence  of 
three-body  and  "switching"  reactions  is  much  like  that  for  NO  in  the  quiet  D-region, 
except  that  the  n  =  1  PH  ion  is  now  generated  first.  Another  important  difference 
is  that  an  intermediate  ion  in  the  sequence,  O^,  reverts  to  upon  reaction 
with  O,  inhibiting  the  overall  ion  transition  scheme  somewhat. 

At  the  top  of  the  D-region,  nominally  90  km,  the  major  positive  ions  are 
usually  O2  and  NO^,  and  negative  ion  concentrations  are  negligible.  Hence, 
the  range  for  the  effective  reconr±)ination  coefficient  is  rather  limited, 

a(02)  <  (90  km)  <  a(NO‘^)  .  (2a) 

The  dissociative  recombination  coefficients  of  these  ions  are  not  very  different. 

For  a  temperature  of  190°  K,  approximate  numerical  values  are: 

3  X  10'^  ''m^  sec"^  <  (90  km.  190°K)  <  6  X  10’^  cm^  sec’^  .  (2b) 

At  a  somewhat  lower  height,  as  discussed  below,  an  ionic  scheme  that  begins  with 
ot  and  NO^  ions  and  ends  with  PH  ions,  H„0^-  (H„0)  becomes  operative.  The 
relative  distribution  of  these  ions  comes  into  thermodynamic  equilibrium  at  some 
lower  height,  roughly  65  km. 

The  transition  altitude,  where  [NO^l  +  [ol]  ==  [HoO^- (H„0)  ,  is  principally 

^  n  w  n 

a  function  of  the  parameters  q.  [  H„0]  and  T.  For  the  quiet  D-region,  this  altitude 

^  12 

generally  is  near  82  km  in  the  daytime  and  near  86  km  at  night  [Narcisi  ].  This 
altitude  becomes  lower,  5-10  km  at  most,  as  q  is  enhanced.  The  transition 
altitude  decreases  as  q  increases  because  with  an  increase  of  q  (and  [e]),  dissocia¬ 
tive  recombination  of  NO^  and  increases,  but  not  their  rate  of  conversion  to 
PH  ions. 


12.  Narcisi,  R.S.  (1974)  Ion  and  neutral  composition  measurements  in  the  lower 
ionosphere  in  Lower  Ionosphere  Studies,  Akademie  Verlag,  Berlin. 
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In  any  case,  the  transition  from  and  o'iNO'*^)  to  aiPII)  is  not  of  very 

much  importance  inasmuch  as  rate  coefficients  for  PH  ions  are  only  4-5  times 
those  for  NO^/O^  ions  and  depend  little  upon  temperature.  Values  of  cVj^  are 

(3.  0  ±  0.  5)  X  10~®  cm^  sec  ^  for  ions  and  about  5  X  lO  ^  cm^  sec  ^ 
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for  the  next  two  heavier  hydrates  [Huang  et  al  ].  Thus,  at  75  km,  where  negative 

“6  3  “1 

ions  are  still  insignificant,  =  4  X  10  cm  sec  is  a  good  estimate  for  most 
conditions.  (The  heights  80  and  85  km  are  near  the  transition  level.  ) 

Negative  ions  become  plentiful  near  70  km,  basically  the  consensus  altitude, 

14 

z(X^),  for  where  A  =  1.  Thomson  scatter  data  (Ganguly  ct  al  ]  agree  with  this 
altitude.  Now  it  is  a  convenient  circumstance  that  the  ion-ion  recombination  coeffi¬ 
cients  of  D-region  ions  neither  differ  substantially  from  one  another  [Smith  et  al^^] 

“*1/2  16 

nor  vary  much  with  temperature,  perhaps  T  [Smith  and  Church^  ],  (Figure  3 
of  the  latter  paper  is  an  excellent  summary  of  and  or.  from  ground  level  to 
90  km. )  Thus,  for  the  lower  half  of  the  D-region,  the  effective  recombination 
coefficient  may  be  expressed  as  (see  Eq.  (1)); 

*  =  (4  X  10'®  +  6  X  10‘®A)(1  +  A)  (3) 

-S3-1  -83-1 

since  Of^  remains  near  4  X  10  cm  sec  and  of^  is  about  6X10  cm'  sec 

This  equation  can  provide  guidance  as  to  A,  if  ^  is  well-determined,  that  is 

A  =  1  where  ^*8X10^  cm^  sec  ^  and  A  =  10  where  4*  *  5  X  10  ^  cm^  sec 

Two  relatively  well-determined  profiles  for  4*  were  deduced  from  the  SPEs  of 

November  1969  and  August  197  2.  The  former  event  was  exceptionally  well  covered 

with  much  in-situ  data  and  both  satellite  and  rocket-based  observations  of  the 

precipitating  proton  fluxes  in  a  program  largely  sponsored  by  DNA.  The  crosses 

in  Figure  1  represent  'I'  for  this  event  from  as  many  as  four  separate  results 

13.  Huang,  C. -M. ,  Whitaker,  M. ,  Biondi,  M.A.,  and  Johnsen,  R.  (1978) 

Electron-temperature  dependence  of  recombination  of  electrons  with 

HgO^- (H20)^- series  ions,  Phys.  Rev.  A.  18:64-67. 

14.  Ganguly,  S. ,  Matthews,  J.D.,  and  Tepley,  C.A.  (1979)  Thomson  scatter  radar 

detection  of  D-region  negative  ions  at  Arecibo,  Geophys.  Rev.  Lett. 

6:89-92 

15.  Smith,  D. ,  Adams,  N.G.,  and  Church,  M.J.  (1976)  Mutual  neutralization 

rates  of  ionospherically  important  ions.  Planet.  Space  Sci.  24:697-703. 

16.  Smith,  D. ,  and  Church,  M.J.  (1977)  Ion-ion  recombination  rates  in  the 

earth's  atmosphere.  Planet  Space  Sci.  25:433-439. 
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[Swider  and  Dean;  Ulwick  ].  Goo<i  agrei,*'nent  with  some  other  results  has  been 
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noted  [Swider  and  Dean  ),  and  the  data  have  been  matched  well  with  a  model 
IQ 

[Swider  et  al  [. 

Circles  in  Figure  1  represent  mean  4*  for  five  instances  dur’ing  the  very  intense 

20 

SPE  of  3-11  August  1972  [Reagan  and  Watt  j.  Electron  concentrations  were 

obtained  with  backscatter  radar.  The  lowest  altitude  data  may  suffer  from  the 

neglect  of  negative  ions  upon  the  radar  signals. 

2  1 

Shownalso  in  the  figure  is  Gledhill's  least  squares  fits  to  all  the  4'  -values 
in  the  literature  for  both  quiet  and  disturbed  (SPE,  SPA,  Flare)  conditions.  These 
two  lines  do  not  differ  substantially  from  one  another  and  underscore  that  the  level 
of  ionization  has  little  bearing  upon  4',  except  to  shift  the  simple  ion/ PH  ion  trans- 
sition  altitude,  which  is  a  function  of  q,  H2O,  and  T.  The  scatter  shown  in  Gled- 
hill's  report,  particularly  in  the  upper  D-region,  is  much  larger  than  can  be 
tolerated  on  the  basis  of  the  chemistry  outlined  above.  In  addition,  a  least  squares 
fit  is  somewhat  misleading  for  the  lower  altitudes  where  a  curved  rather  than  a 
linear  fit  is  undoubtedly  proper.  Attachment  becomes  more  important  with  decreas¬ 
ing  altitude  and  the  initial  process  is  dependent  on  the  square  of  the  neutral  con¬ 
centration.  The  divergence  of  these  curves  below  65  km  is  questionable.  On  the 
other  hand,  the  higher  4*  values  above  70-75  km  for  the  quiet  curve  is  consistent 
with  the  notion  that  4*  increases  as  q  decreases  since  dissociative  recombination  of 
electrons  with  PH  ions  rather  than  or  NO^  is  more  probable  as  these  ions  have 
a  greater  chance  to  transition  to  PH  ions  prior  to  dissociative  recombination  when 
[e]  is  low. 
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Figure  1.  Effective  Recombination  Coefficients  for  the  November  1969, 
SPE  (crosses)  [Swider  and  Dean)**^  and  the  August  1972,  SPE  [Reagan  and 
Watt^O]  (circles)  Compared  to  a  Least  Squares  Fit  of  All  Such  Data  for 
Quiet  and  Disturbed  Conditions  [Gledhill^l] 


-  R  -.1 

Recalling  from  Eq.  (3)  that  'I' -  8  X  10  cm  sec  at  z(A  j',  the  four  curves 

in  Figure  1  yield  a  range  of  about  63  to  68  km  for  z(A  j).  It  seems  reasonable  to 

assume  that  this  altitude  should  decrease  somewhat  as  q  (and  [e])  increases,  since 

the  probability  of  an  electron  undergoing  dissociative  recombination  must  increase 

vs  attachment  which  is  basically  independent  of  fe).  However,  the  change  of 

19 

z(A  ^)  is  not  great.  Recall  that  Swider  et  al  were  able  to  model  the  November  1969, 

SPE  rather  well.  Their  z(A  .)  was  62  km.  The  same  model  yielded  63  km  for 

^  22 

z(A  j)  under  quiet  conditions  [Swider  ).  However,  both  altitudes  appear  to  be  a 
little  too  low. 


22,  Swider,  W.  (1983)  Latitudinal  influences  on  the  quiet  daytime  D-region, 
Adv.  Space  Res.  2:213-216. 
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Profiles  of  the  effective  recombination  coefficient  determined  for  quiet  daytime 

22 

conditions  using  several  different  atmospheric  (T,  M)  models  [Swider  )  are 
compared  in  Figure  2  to  the  SPE  data  shown  in  Figure  1.  A  larger  ♦  for  quiet 
conditions  as  compared  to  the  November  1969,  SPE  is  expected  above  about  7  5  km 
as  already  stated.  It  is  somewhat  surprising  that  'J' from  SPE'72  fits  the  quiet 
(low  q)  model  so  well.  These  results  emphasize  that  ^  is  a  function  of  [M],  the 
total  concentration.  This  is  the  reason  for  the  shift  among  the  curves  with  altitude. 
There  is  no  difference  between  curves  1  and  4  (same  atmospheric  model)  except 
above  70  km,  where  their  NO  distributions  are  different.  In  any  case,  the  two 
SPE  data  sets  appear  to  encompass  both  disturbed  and  quiet  conditions. 


1- 1976  U.  S.  STAND.  "NORMAL"  lxlO‘'^M] 

2- CIRA'72-60*N(N0V)  ••  2xI0‘'1m] 

3- CIRA'72  -  70®N(JUL)  "  2xI0‘'^[M] 

4-  1976  U  S.  STAND.  LOW  2xI0''®[M] 

HoO  -  ALLEN  et  al.  (1981)-  MODEL  3 


EFFECTIVE  RECOMBINATION  COEFFICIENT  <)^=q/[e]  (cm’s"') 


Figure  2.  Theoretical  Profiles  of  Effective  Recombination  Coefficients 
[Swider22]  for  Several  Atmospheric  Models  as  Compared  to  Two  SPE 
Results  From  the  Previous  Figure 


An  example  of  the  use  of  an  empirical  •i  may  be  illustrated  for  the  23  September 
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197  8,  SPE.  Using  the  q -values  listed  by  Kossey  et  al  for  noon,  and  the  'i  values 

17  4  -3 

of  Swider  and  Dean,  the  following  [e]  are  deduced:  2,  9  X  10  cm  (80  km), 

1.  4  X  10^  cm^  (70  km),  5.  1  X  10*'  cm  ^  (60  km),  4.  0  X  10^  cm^  (50  km)  and  about 
-3 

2.0  cm  at  40  km.  These  concentrations  were  used  to  calculate  the  electron 
conductivity  from  the  formula 


a 

e 


3.  27  X  lO'^'^te] 
- P^Tor'r) - 


mho/m  . 


(4) 


The  results  are  plotted  (dots)  in  Figure  3  along  with  the  results  obtained  from  the 

(extrapolated)  inversion  of  the  data  and  from  a  detailed  model.  Clearly,  this 

simple  method  of  determining  a„  from  an  experimentally  derived  'J'  which  is 
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backed  by  physical  theory  [Swider  et  al  ].  works  well. 


CONDUCTIVITY  (mho/m) 


Figure  3.  Conductivity  as 
Determined  from  a 
Detailed  Model  and  t  ie 
Inversion  of  Propag.  ion 
Data  [Kossey  et  al],  as 
Compared  to  a  Simpb' 
Estimate  (dots)  Derived 
From  Effective  Recombination 
Coefficients 


23.  Kossey,  P.  A.  ,  Rasmussen,  J.  E.  ,  Field,  E.  C. ,  and  Warber,  C.R.  (1985) 

Conductivity  profiles  of  the  disturbed  polar  ionosphere  from  VLF  reflection 
data,  AGARD-CP-382:8.  1.  1-8.  1.  10. 
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3.  EFFECTIVE  RECOMBINATION  COEFnCIENT  AT  NIGHT 

Above  80  km,  there  is  little  difference  between  night  and  day  for  disturbed 

conditions  (same  q  day  and  night).  For  quiet  conditions,  q(day)  >>  q(night)  and 

the  transition  altitude  from  simple  diatomic  ions  to  PH  ions  is  higher  at  night 
12 

[Narcisi  ].  Negative  ion  formation  is  minimal  both  day  and  night  above  80  km 
where  there  is  considerable  [O]  to  counter  the  formation  of  the  initial  negative  ion, 
O^,  by  the  process  O  +  ©2  "*■03  ^ 

At  night,  photo-processes  are  inopera*ive  and  atomic  oxygen  declines  sharply 
with  decreasing  height  below  80  km.  Attachment  rapidly  becomes  the  main 
electron  loss  process  with  decreasing  altitude  below  this  point.  There  is  a  quick 

switch  from  ^  ol  (z  >  80  km)  to  'I<  =  Xo.d  +  X  )  for  z  <  70  km.  It  is  convenient 

^  '  -30 

to  express  attachment  in  terms  of  the  0„  formation  rate,  L(A)  »  1.8X10 


[0„]^  sec“\ 


The  formula  [Swider  et  al^^] 
V  +  -  ( L(A) 


appears  to  be  quite  adequate  for  disturbed  conditions.  The  equation  simplifies  to 

[e]  =  (q/Oj^)  near/above  80  km,  that  aS,  <!»  -  Q^,  and[el  =  q/L(A)  near/below  70  km 

vihich  yields  4*=  L(A)^/q.  Thus,  at  night,  a  pure  attachment  (linear)  loss  of 

electrons  occurs  below  this  point.  The  electron  concentration  is  proportional  to 

the  production  rate.  Gledhill  calculated  a  mean  least  squares  fit  for  4*  under 

nighttime  conditions.  This  result  should  not  be  used  in  view  of  Eq.  (5).  The  quiet 

nighttime  D-region  is  basically  nonexistent,  being  only  the  tail  of  the  lower  night- 
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time  E -region  [Keneshea  ].  However,  Eq.  (5)  should  be  appropriate  for  this 

-7  3  - 1 

case  too.  The  recommended  value  for  aj^  is  6  X  10  cm  sec 


24,  Swider,  W.  ,  Narcisi,  R.  S.  ,  Keneshea,  T.  J. ,  and  Ulwick,  J.  C.  (197  1) 
Electron  loss  during  a  nighttime  PCA  event,  J.  Geophys.  Res.  76:4691-4694. 

25.  Keneshea,  T.  J,  ,  Narcisi,  R.  S.  ,  and  Swider,  W.  (1970)  Diurnal  model  of  the 

D-region,  J.  Geophys.  Res.  75:845-854. 
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4.  DISCUSSION  AND  CONCLUSIONS 


Distributions  of  [O],  [NO],  [H2O],  and  T  are  important  factors  in  formulating 
physical  models  of  the  D-region.  Global  models  of  these  parameters  are  mostly 
incomplete.  Other  species  distributions  relevant  to  the  negative  ion  chemistry,  like 

NO„,  and  0„,  are  needed  also.  Rate  constants  for  many  of  the  positive  and  nega- 

^  26 
tive  ion  reactions  are  unmeasured.  There  is  evidence  [Narcisi  et  al  ]  that  a 

27 

layer  of  heavy  negative  ions  exists  near  80  km.  Backscatter  data  [Ganguly  ] 
yield  similar  evidence  about  a  thira  of  the  time.  Gas -phase  chemistry  allows 
for  only  a  few  light  negative  ions  at  this  altitude. 

These  and  other  complexities  not  withstanding,  the  D-region  is  generally 
understood.  Physical  models,  while  imperfect,  provide  the  basis  of  this  under¬ 
standing.  However,  to  provide  a  predictive  model  of  electron  concentrations  that 
can  match  propagation  data,  parameters  like  NO,  H2O,  and  T  need  be  observed 
virtually  in  real  time.  Constant  monitoring  of  these  parameters  is  impractical, 
although  real  time  monitoring  of  the  particle/X-ray  intensity  is  required  for 
disturbed  events.  On  the  other  hand,  the  effective  recombination  coefficient  for 
electrons  in  the  daytime  D-region  appears  to  be  much  less  variable  than  either 
[e]  or  q. 

2 1 

The  considerable  scatter  for  ^  in  the  literature  [Gledhill  ]  probably  results 

from  poor  measurements  of  either  [ej,  q,  or  both  as  compared  to  better  measured 

events  like  the  SPE's  of  November  1969  and  August  1972,  The  parameter  4*  cannot 

vary  that  widely  on  physical  grounds.  At  65-85  km  in  the  daytime,  electrons 

dissociatively  recombine  with  either  simple  ions  (NO^,  ot)  or  PH  ions, 

+  ^ 

HgO  •  (H2O)  .  The  latter  ions  become  proportionally  more  important  with  descent. 

A  knowledge  of  T  and  H2O  can  help  refine  the  value  of  although  [e]  only  depends 

upon  the  square  root  of  'I'  »  0!|^  anyway  at  these  altitudes.  The  main  effect  of  a 

larger  q  is  to  lower  the  crossover  altitude  of  simple  ions  to  PH  ions  a  few  kilometers 

since  with  a  greater  number  of  electrons,  dissociative  recombination  with  Ot  and 
+ 

NO  ions  is  enhanced  relative  to  their  transition  to  oxonium  ions  which  have  larger 
dissociative  recombination  coefficients.  At  lower  altitudes,  where  negative  ions 
are  significant,  ’J<  is  perhaps  more  nearly  invariant  (Figures  1,  2,  and  3). 


26.  Narcisi,  R.,  Bailey,  A.,  Federico,  G. ,  Wlodyka,  L.  (1983)  Positive  and 

negative  ion  composition  measurements  in  the  D-  and  E -region  during 
the  26  February  1979  solar  eclipse,  J.  Atmos.  Terr.  Phys.  45:461-478. 

27.  Ganguly,  S.  (1984)  Observation  of  large  negative  ion  zones  around  88  km, 

J.  Atmos.  Terr.  Phys.  46:633-636. 
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Additional  data  are  needed  to  test  this  hypothesis.  A  DNA-sponsored  program  is 

underway  to  do  just  that.  Electrons  are  quite  difficult  to  measure  at  low  altitudes, 

but  a  new  technique  will  be  tried.  Ambient  electrons  will  be  converted  to  negative 

ions  via  a  release  of  SF„.  Non-ambient  SF.  ions  will  be  measured  by  a  mass 

o  t> 

spectrometer  relative  to  the  ambient  negative  ion  species. 

A  daytime  profile  of  'i'  for  50  to  90  km  averaged  from  the  crosses  and  circles 
shown  in  Figure  1  is  recommended.  The  circles  may  be  preferable  for  lower  T 
(summer)  and  the  crosses  for  winter.  The  circles  also  appear  to  represent  quiet 
conditions  best  (see  Figure  2).  At  night,  is  not  invariant,  being  inversely  propor 
tional  to  q  below  about  70  km.  However,  a  suitable  analytical  solution  exists  for 
[e]  under  disturbed  conditions,  which  also  should  be  applicable  for  quiet  conditions. 


^  WWWWfcfVif. 
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